Spectral correlations: understanding oscillatory contributions 
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We give a transparent derivation of a relation obtained using a supersymmetric non-linear sigma 
model by Andreev and Altshuler [Phys. Rev. Lett. 72, 902, (1995)], which connects smooth and 
oscillatory components of spectral correlation functions. We show that their result is not specific to 
random matrix theory. Also, we show that despite an apparent contradiction, the results obtained 
using their formula are consistent with earlier perspectives on random matrix models. In particular, 
the concept of resurgence is not required. 



Spectral correlations of complex quantum systems, 
such as disordered metals and classically chaotic quan- 
tum systems, are known to be nearly universal. For small 
ranges of energy they are well approximated by the Gaus- 
sian invariant ensembles of random matrix theory (GXE, 
where X=0,U or S stands for orthogonal, unitary and 
symplectic invariance) Deviations from GXE be- 

haviour at larger energy scales may be consistently incor- 

Pted using semiclassical or perturbative approaches 
. An interesting development in this field was a pa- 
per by Andreev and Altshuler (AA) ||], who introduced 
a relation which suggests a degree of non-universality in 
short ranged spectral correlations. Their calculations are 
based on the non-linear sigma model. Our paper will give 
a transparent physical insight into their relation. 

AA 1^ consider the spectral two-point correlation func- 
tion, defined as 



Rf){e) = {d{E + e/2) d{E - e/2)) - 1 



(1) 



Here d{E) = J2n^i-^ ^ -^") ^^e density of states, i?„ 
are the eigenvalues of a Hamiltonian H and A(i?) — 
{d{E))~^ is the mean level spacing which is assumed to be 
independent of E in the following. The index f3 = 1,2,4 
distinguishes orthogonal, unitary or symplectic symme- 
try classes respectively. 

AA divide Eq. (|l]) into a smooth and an oscillatory 
contribution, and propose that (for (3 — 2) these are re- 
lated as follows 
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2 cos(27re/A) £>(£), 



(2) 
(3) 



The relation is exact for the GUE with a suitable choice 
for -D(e) and generalizations are approximately true for 
the other ensembles. AA propose that Eqs. (2) and (3) 
are quite general, and should be used to predict the os- 
cillatory component from the average component. The 
quantity D{e) is a spectral determinant. 

The AA relation contained in (||) and (||) has attracted 
considerable attention (see, for example, @-||]), partly 



because of the suggestion that it contains information 
about 'resurgence'. This is to be interpreted in terms of 
Gutzwiller's [ p^ relationship between periodic orbits and 
the density of states: periodic orbits orbits with period 
tj are associated with oscillations in the density of states 
of period ej ~ Ti/tj. The concept of resurgence is that 
information about long period orbits is encoded in prop- 
erties of the short orbits. In this context, W"{e) may 
be derived using Gutzwiller's relation from properties of 
short orbits, and the AA relation then gives information 
about fiuctuations in the spectrum with wavelength equal 
to the mean level spacing A, corresponding to orbits of 
period equal to the 'Heisenberg time', t-a = 27r?i/A. 

This interpretation of the AA result is at variance with 
some old ideas about the applicability of random matrix 
theory. Dyson introduced a Brownian motion model 
for random matrix theory, giving a Langevin equation of 
motion for the response of energy levels to a stochastic 
perturbation. Dyson later suggested [|l2| that this model 
may give valuable insights into why random matrix the- 
ory applies to generic systems. He considered the dy- 
namics of the energy levels under the effect of a stochas- 
tic perturbation, which is large enough to shift energy 
levels significantly, but which remains small enough to 
leave other features of the system (such as the classical 
dynamics) unchanged. Dyson showed that the equations 
of motion simplify if the discrete Fourier coefficients of 
the level displacements are used as dynamical variables: 
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^ AEn exp{-27rikn/N) 



(4) 



n=0 



where AEn — En~nA, N is the number of energy levels, 
and k takes N successive integer values. We will take the 
maximal k to be int(A^/2). The long wavelength modes 
evolve almost independently, with a long relaxation time, 
which scales as k~^ |12|1]- The short wavelength modes 
remain strongly coupled. The stochastic perturbation 
will bring the short wavelength modes into equilibrium, 
giving statistics of the Uk which are identical to the ran- 
dom matrix ensemble. 
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These arguments were later supported and extended 
with the aid of semiclassical estimates of matrix elements 
Q. The argument in assumes that the matrix ele- 
ments of the perturbation are uncorrelated. It was shown 
that semiclassical estimates are consistent with this hy- 
pothesis when considering the stochastic force driving the 
large k modes. The forces driving the small k modes are 
modified by the classical dynamics of the system. The 
resulting picture is that long wavelength fluctuations are 
non-universal, but that at short wavelengths the exci- 
tations of the modes are precisely the same as for the 
appropriate GXE. In particular, there is no modification 
of the statistics of the Fourier coefficients at unless \k\/N 
is small. This appears to be in contradiction with the AA 
relations, which suggest that non-universal corrections to 
the smooth part of the correlation function are echoed by 
oscillations with large wave numbers k. 

One resolution of this contradiction would be that 
there are previously un-suspected correlations between 
matrix elements due to 'resurgence', which are not cap- 
tured by the semiclassical approximations in |^|. We will 
however show that there is no need in invoke this princi- 
ple. We will first describe a simple derivation of the AA 
result, and comment on its applicability. We will then 
describe the Brownian motion model and use it to derive 
the correlation function Rp{e) for the case of a system 
with diffusive electron motion. In this calculation we 
assume that the statistics of the Fourier coefficients are 
unchanged for large \k\. The fact that we reproduce ex- 
isting results verifies that resurgence is not an essential 
ingredient in their derivation. Finally we will comment 
on the correlation function for classically chaotic systems. 

Our starting point is the following general expression 
for the correlation function i?^(e): 



upon expanding the slowly varying window function in 



n— — OO 



(5) 



where p^(7i, e)de denotes the probability of finding that 
the difference between Eq and £"„ is in the interval 
[e, e -|- de]. We will show below that in many cases, the 
P/3{n, e) are well approximated by Gaussians for large n. 
We will estimate the variance <Jp{n) for different systems. 

First, however, we show how the relations (||) and (j^) 
can be derived from (S). We write 



Here i?«^(e) is defined as 



de'w{e-e')Rp{e') 



(6) 



(7) 



where ^(e) is a suitable window function (which could 
be a Gaussian centred around zero with variance much 
larger than ap{n) and normalized to A^^). R^^'^{e) is the 
remaining oscillatory contribution. Using (ph we have. 



(0) 



R^^e) = -1 + AY, / de'w{e~e')p0{n,e') (8) 

oo ^ 

~ -1 + A ^ [w{e-nA) + -w" {e-nA)al{n)] . 
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Assuming that cr|(ri) is a slowly varying function of n, 
we may approximate the sum over n as an integral and 
obtain 



2^1^ -^H, ^-nA. (9) 



Consider now the oscillatory contribution R°^'^{e). Again 
assuming that cr^ (n) is slowly varying, these can be eval- 
uated from Eqs. (|^) using Poisson summation: 

Rf^{e) 2cos(27rMn) e-^2n^,)''^l{n)/2A^ ^^^^ 

This sum is dominated by the fj. = 1 term. Defining 
D{e) = 47r^ exp[— 27r^o-^(n)/A^] we obtain the relations 
(^ and (^ . We remark that these formulae are valid for 
any correlation function that can be approximated as a 
sum of Gaussians, in regions where the variance satisfies 
cr^ ^ A^, and where varies sufficiently slowly as a 
function of n. Relating D{e) to (T^(n) gives a clear in- 
sight into its meaning. The relations (g) and (||) may be 
easily extended by considering higher-order derivatives in 
(^ and higher Fourier components in (^). Extensions to 
non Gaussian spacing distributions are possible. Also, we 
emphasize that these relations are not specific to spectral 
correlation functions. For example, they are applicable 
to density correlations in solids: for thermal excitation of 
phonons the corresponding two-point function is a sum 
of Gaussians, and a^{n) ~ logn in d = 2 dimensions p3[ . 

Next we consider the calculation oi pi3{n, e) and crj(n) 
using Dyson's Brownian motion model. In the Brownian 
motion model, the matrix elements of the Hamiltonian 
H undergo a diffusive evolution as a function of a ficti- 
tious time variable r. We denote the infinitesimal change 
of H by 5H. In the GOE case, we have for n > m and 
n' > m' 



(SH,nn) — , {5Hmn SH„ 



(11) 



(extensions for GUE and GSE are given in ||ri|). The 
diagonal elements obey 

{SH„n)=0, {6H^raSH„n) =2STr^C^^S. (12) 

Dyson |l^ originally discussed the case where G^nn'^ = 
Smn and = 1, for which the statistics of the En 

are the same as for the GXE. We will argue below that 
non-universal deviations from the GXE are encoded in 
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C^'!*;, (and C^^ = C°„"_„). Using second order 



perturbation theory leads to a Langevin equation 



= J2 



5H„ 



for the energy level shifts 5 En- Thus 



{SE,ndEn)=2STr^C^t^ 



(13) 

(14) 
(15) 



Semiclassical estimates indicate that the de- 
crease for large values of n, i.e. the repulsive interaction 
is screened at long range. This effect was considered in 
|l5[ ; for our purposes it is not significant, and we will set 
= 1 throughout. We now use (jj) to obtain the same 
equations of motion in terms of the Fourier variables at- 
Using SAE„ = 5 En, we have 



(16) 



where h = iV"^ C'^''''^ exp(-27rifcn/iV). The trans- 
formation of the drift term is less straightforward. In 
general the expectation value of 5ak is a complicated 
function of all of the Ofc, but in the limit \k\/N 
the equations decouple and obey [| 11 



(17) 



By solving a Fokker-Planck equation for the real and 
imaginary parts of a^, these are found to have a steady 
state distribution which is Gaussian. The value of (| a^p) 
can be deduced by requiring that (|afc(T -I- ^t)P) = 
(|afc(T)p): writing afc(T -|- 5t) — akir) + Sak and using 
(P) and (|l^) we find, provided |fc|/iV < 1: 



\ak\ 



2TT^I3k 



(18) 



We cannot determine (|afcP) for larger values of |fc| from 
this approach. We now make our key assumption, that 
Ik ^ N^^ when \k\/N is not small. This corresponds 
to assuming that there are no short-ranged correlations 
between the diagonal matrix elements, i.e. C^'J^^ ~ ^nm- 
Semiclassical arguments which support, but do not prove, 
this assumption are given in ^ . (The existence of short 
ranged correlations would represent a type of 'resur- 
gence', in the sense discussed earlier). Under this as- 
sumption, the equations of motion of the Uk are identical 
to those for the Brownian motion model describing Gaus- 
sian invariant ensembles. We therefore conclude that 
when |fc|/7V is not small, the mode intensities (|afcp) are 
identical to those of the Gaussian invariant ensembles. 

For large n, the level spacings £"„ — Eg are seen to be a 
sum of many independent random variables and are thus 
Gaussian distributed 



with mean nA and with variance 



sin^(^) 



(20) 



Writing t = t^k/N, and Nik = I{t), Eq. (|2g) becomes 
4A2 



f3n 



2 /o T''''^')^^"^^ 



(21) 



where Ji3{t) ~ I{t) for t ^ 1, and Jpit) takes a univer- 
sal (but to us, unknown) form when t is not small. For 
the Gaussian invariant ensembles (where I{t) = 1), the 
variances clearly grow logarithmically with ri. They are, 
asymptotically for large n, 



2A2 



[log(2^n) + C0] 



(22) 



with C'l = -log\/2, C2 = and C4 = log(4/7r). Using 
this expression in (^ and (0) gives the correct leading 
order contributions to R^(e) and R'^^'^{e) in the limit 
e/A^oo (c.f. @). 

Next we consider how the function I{t) must be modi- 
fied at small t to take account of classical dynamics. For 
small values of the amplitude Sak can be estimated 
semiclassically [0 



Sak ~ ^tv[6HU{t)] 



t = 



2Trh k 
A~ N 



k 
N 



(23) 



where U{t) — exp{—iHt) is the evolution operator. We 
consider first diffusive systems (electrons in disordered 
metals), then systems with a chaotic classical limit. 

Diffusive systems. In this case we may consider the 
perturbation SH to be uncorrelated random changes of 
the site energies Un in an Anderson tight-binding model 



5iJ = ^(5U„Pn, (SV^SV^,) 



(24) 



where Pn is the projection for locating an electron on lat- 
tice site n. Using the semiclassical approximation ([2^), 
I{t) is seen to to be proportional to the probability of 
returning to the original site after time t. Normalising so 
that I{t) approaches unity for large t, we have 



-Dktt 



(25) 



where the sum is over the eigenmodes of the Helmholtz 
equation (V^ -I- k1)^^{v) = with Neumann boundary 
conditions. In a quasi-one dimensional system, k^ = 
TxvjL. In this case we obtain from Eq. (Elh 
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al{n) 



2A2 



/37r2 



log(27rn) + C^ + -^log(l- 



(26) 



where g = TT^fiD/ L'^A is a dimensionless conductance. 

An alternative way to express the results for Rpie) is in 
terms of the 'form factor' K{t) which is the Fourier trans- 
form of Rpit)- As discussed by AA, Eq. (^ gives rise to 
non-universal structures in the form factor at the Heisen- 
berg time in- We have verified their existence in numer- 
ical simulations: we used an ensemble of complex Her- 
mitian banded random matrices (of dimension N = 1000 
and band width h = 35) modelling a quasi-one dimen- 
sional diffusive system. We fitted the dimensionless con- 
ductance, using states from the centre of the spectrum 
(obtaining g ~ 2.0). The results are shown in Fig. |l| and 
are in good agreement with the theoretical predictions, 
Eqs. (9), (10) and (26). 

The theoretical predictions shown in Fig. 1 are pre- 
cisely equivalent to the AA results (7) and (14) in [||. 
We have thus shown that their results are consistent with 
the much simpler approach to justifying random matrix 
theory discussed in p^ . 

Classically chaotic systems. Consider the expression 
( ^ ) for the fluctuation 5ak of the Fourier coefficients. 
If the system has a smooth classical Hamiltonian, this 
expression will clearly be negligible unless t corresponds 
approximately with the period of a periodic orbit. More- 
over, it was shown in jij that when the motion is chaotic 
the 5ak have statistics corresponding to random matrix 
theory for large t. A simple model capturing the essential 
features of this case is to replace the lower limit of the 
integral in ( |2l| ) with the period of the shortest periodic 
orbit. This gives: 



2A2 



a(?^)-,„,(i.)-,.c. 



(27) 



This model has the feature that tT| (n) is finite in the limit 
n oo, corresponding to the behaviour of Dyson's A3 
statistic for systems with a smooth classical limit |3| . We 
note that this implies [using (||) and (^] that the oscil- 
latory part of the correlation function does not decay to 
zero. The form factor is seen to have a delta function at 
the Heisenberg time with a magnitude proportional 
to (to/^n)^^''- This feature has not been remarked upon 
in earlier papers which have discussed the form factor for 
chaotic systems 

A more precise estimate of u^in) for specific systems 



can be obtained usin 
approach used in 



periodic orbit theory, following the 
The conclusions are unchanged: 
(n) remains finite as n — > cxd, and there must exist 
oscillations in the correlation function which do not de- 
cay. One difference is that for large n the pp{n,e) are 
dominated by the shortest orbits. Because only a finite 
number of components are significant, the central limit 
theorem cannot be used to assert that these distributions 
are Gaussian. Moreover small deviations from a Gaus- 
sian distribution can have a large effect on the Fourier 



transform of the distribution, which determines the mag- 
nitude of the oscillations in (|lo|). We infer that the AA 
relations may need to be modified when applied to the 
spectra of systems with a smooth classical limit. 

Conclusions. We have shown that the AA relations 
have a simple interpretation, independent of random ma- 
trix theory and supersymmetric techniques. We have 
shown how spectral correlations of diffusive systems are 
obtained using from Dyson's Brownian motion model, 
indicating that the AA relation is consistent with this 
approach. 
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FIG. 1. The form factor K{t) for a quasi-one dimensional 
diffusive system. Numerical simulations for banded random 
matrices (•), are compared with theoretical results ( ) ac- 
cording to Eqs. (|), (|l|) and (||). The CUE result is also 
shown ( ). 
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